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Through s — d coupling, a superconducting thin film interfaced to a noncollinear magnetic insulator 
inherits its magnetic order, which may induce unconventional superconductivity that hosts Majorana 
edge states. From the cycloidal, helical, or (tilted) conical magnetic order of multiferroics, or the 
Bloch and Neel domain walls of ferromagnetic insulators, the induced pairing is {px + P 2 /)"Wave, 
a pairing state that supports Majorana edge modes without adjusting the chemical potential. In 
this setup, the Majorana states can be separated over the distance of the long range magnetic 
order, which may reach macroscopic scale. A skyrmion spin texture, on the other hand, induces a 
(jpr + zp(^)-wave-like state, which albeit nonuniform and influenced by an emergent electromagnetic 
field, hosts both a bulk persistent current and a topological edge current. 
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Introduction.- Motivated by possible applications in 
non-abelian quantum computation [T] , the search for Ma¬ 
jorana fermions in condensed matter systems has wit¬ 
nessed a boost recently lai]. Indeed there has been 
much effort to design and fabricate one-dimensional het¬ 
erostructures in which topological p-wave superconduc¬ 
tivity is proximity induced HHig. One particularly 
promising proposal are chains of magnetic atoms with 
noncollinear spin texture on the surface of a conventional 
superconductor (SC) pMll22] . The presence of these mag¬ 
netic adatoms induces Shiba bound states [23H2S] , whose 
low-energy physics is equivalent to a one-dimensional 
(ID) p-wave SC with Majorana modes at its ends. Scan¬ 
ning tunneling measurements of zero-bias peaks at the 
ends of Fe chains deposited on superconducting Pb have 
been interpreted as evidence of Majorana modes [26], al¬ 
though no general consensus has been reached regarding 
the definitive existence of Majorana states in these sys¬ 
tems [27]. 

Since thin films are generally easier to manufacture 
than adatoms, it is intriguing to ask if such ID pro¬ 
posals can be generalized to two-dimensional (2D) sys¬ 
tems, where a superconducting thin film is coupled to a 
noncollinear magnet. The work by Nakosai et al. [28] 
first shed light on this issue. It was found that a spiral 
spin texture in proximity to an 5-wave superconductor 
induces a {px +Py)-wave state, while a skyrmion crystal 
spin configuration gives rise to {px + ip^)-wave-like pair¬ 
ing. The former paring state exhibits bulk nodes and 
Majorana fiat-band edge states, whereas the latter one is 
characterized by a full gap with chirally dispersing Ma¬ 
jorana edge states that carry a quantized Hall current. 
In this Letter, we show that these phenomena occur in a 
much broader class of superconductor/noncollinear mag¬ 
net interfaces. In particular, we find that a great vari¬ 
ety of noncollinear magnets, including multiferroic insu¬ 
lators with helical, cycloidal, and (tilted) conical order, 
as well as magnetic domain walls of ferromagnetic (FM) 


insulators, interfaced with an 5-wave superconductor in¬ 
duce a {px + Py)-wave pairing state with Majorana fiat 
bands at the boundary. We derive the general criterion 
for the Majorana edge state for any wave length and di¬ 
rection of the noncollinear magnetic order, show that the 
Majorana modes occur without fine-tuning of the chem¬ 
ical potential, and demonstrate that the translation in¬ 
variance along the direction perpendicular to the non¬ 
collinear order greatly enhances the chance to observe 
Majorana states. Furthermore, we investigate a single 
skyrmion spin texture coupled to an 5-wave superconduc¬ 
tor and shown that at this interface an inhomogeneous 
{pr + ipc^)-wave-like pairing is induced, which coexists 
with the emerging electromagnetic field resulted from the 
noncoplanar spin texture. 

SC/multiferroic interface.- Evidence from the Tc re¬ 
duction in magnetic oxide/SC heterostructures [29] [30] 
suggests that s — d coupling FS^ • a generally exists at 
the interface atomic layer between an SC and an insu¬ 
lating magnetic oxide [31], where cr is the conduction 
electron spin and the local moment. This leads us 
to consider the following model for the SC/multiferroic 
interface, which is the 2D generalization of the ID pro¬ 
posals of Refs. Iini [HI 


H = 


i,S,a 


fiafi-\-6a ^iS fd' 


/ V la 


fiJu 


+ E • <y)a0 flfip + E ^0 fyi + FP) > 
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where i = {ix/y} is the site index, S = |a, b| is the 
planar unit vector, and a is the spin index. In the ab¬ 
sence of spin-orbit interaction, the majority of the non¬ 
collinear order discovered in multiferroics, for instance in 
perovskite rare earth manganites [32H35] . can be generi- 
cally described by the conical order 

Bi = = (B|| cosOi) , (2) 
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as far as their effect on the SC is concerned, since the 
choice of coordinate for cr is arbitrary. For instance, cy¬ 
cloidal and helical order are equivalent by trivially ex¬ 
change two components of cr, and conical orders with 
any tilting angle are equivalent. Here B\\ = r|S^||| and 
= rjSi^l are the planar and out-of-plane components 
of the local moment, and the planar angle 6 ^^ = Q • 
is determined by the spiral wave vector Q and the pla¬ 
nar position r^. The cycloidal order is the case when 
B±_y = 0. Note that a three-dimensional conical order 
projected to a surface cleaved at any direction is still that 
described by Eq. Q, so our formalism is applicable to a 
thin film or the surface of a single crystal multiferroic in 
any crystalline orientation, assuming no lattice mismatch 
with the SC and a constant |S^| = |S|. 

We perform two consecutive rotations to align the 
field along cr^, 

( fit \ ^ ( cosf -sinf \ ( cos| zsin| \ f gn\ 

\ fa J V sin f- cos f- J \ ism^ cosl J \ ) 


where siny = B^/Bq and Bq = + B\ = |rS|, 

under which one obtains 
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In what follows, we consider hopping to be real and 
isotropic tis = t'^s ~ 

In the limit Bq ^ |/i| ^ {t,Ao} with /i < 0 , one 
can construct an effective low energy theory for the spin 
species near the Fermi level, which is the spin down band. 
This is done by introducing a unitary transformation 
H' = to eliminate the spin mixing terms order 

by order [13]. At first order, the pairing part 


^eff,A = 

i,S 



^ ^otPisgn.gi+si + h.c. 


( 5 ) 


resembles a spinless p-wave superconductor with 
anisotropic nearest-neighbor and next-nearest-neighbor 
hopping. From Eq. Q one has = cos 7 sin Q • ^/2, so 
the induced pairing is of {px +Py)-wave symmetry, with 
the magnitude of the gap determined by the wave length 
of the planar component of the conical order. 


To derive the criterion for the appearance of Ma- 
jorana edge states, we introduce Majorana fermions 
gia = +ibi 2 cr), gL = ^ {biia - ibi 2 a) with 

{bima,bi'm'0} = 25u'5mm'ba0, and express the Hamil- 
tonian B[ = (V4) ^q^(q)b-q in terms of the basis 

= (^qlt5^qU5^q2t) &q24,)[Ill- We choose open bound¬ 
ary condition (OBC) along x and periodic one (PBC) 
along y, such that —tt < Qy < tt is a good quantum 
number. At a particular qy, only when satisfies 

(3a sin qx + Pb sin % = 0 , ( 6 ) 

is the Majorana Hamiltonian A^q^^qy) skew symmetric. 
The two solutions q^i and qx 2 are the high symmetry 
points at which the Pfaffian is calculated. The topologi¬ 
cal index [14] is now a function of qy 

= Sign (Pf[^(fei,gj^)]) Sign (Pf[^(fe 2 ,gy)]) • (7) 

When M{qy) = —1, or equivalently 

\l^1 + (Im - 2tabCosqy \ + \2taaCosqxi\f > B 

> \l ^1 + {\g - 2tabCosqy\ - \2taa cos qxi\f , ( 8 ) 

where as = {as + /2, the Majorana edge state with 

momentum qy appears. Setting as = 0 and qxi = 0 
recovers the well known ID result [T4] . 

Equation is the general criterion for the Majorana 
state to appear at momentum qy for any given spiral or 
conical order. It is supported by numerically solving the 
Bogoliubov-de Gennes (BdG) equation with the bound¬ 
ary conditions we choose. A spin-generalized Bogoliubov 
transformation gi^ = Y.n^incr»lna+v*naalla is intro¬ 
duced to diagonalize Eq. Q. Figure [^b) shows a typical 
dispersion E{n,qy), which display Majorana zero-energy 
states in the g^’s that satisfy Eq. Note that Equation 
^ can be satisfied even if /i = 0 , so adjusting chemical 
potential is generally not needed. Consequently, the edge 
states can occur in an isolated sample without attaching 
any leads. The 2tas cos qy factor greatly enlarges the 
number of g^’s that can satisfy Eq. ([^, hence increases 
the chance to observe Majorana Fermions, as one can 
see from the phase diagram shown in Fig. [^a) that has 
much larger topologically nontrivial phase in the Bo-a 
space than the ID cases, where the (weak) topological 
phase [36] is judged by whether any qy satisfies Eq. ^ 
at a given (Hq, a). 

The localized Majorana edge states can be seen as 
the zero bias peaks (ZBPs) in the local density of states 
(LDOS) along x direction, as plotted in Fig. [^c). Three 
gap-like features show up in the LDOS. The two sym¬ 
metric in uj come from the bulk gap A whose position is 
shifted by the s — d coupling ±Bo as its effect is similar 
to a magnetic field, and typically has a magnitude [37] 
0.01 ^ leV so Ho > A. The gap-like feature near zero en¬ 
ergy represents the induced (pa,-hpy)-wave gap in Eq. 
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FIG. 1: (color online) (a) The topologically nontrivial phase 
(blue region) of the 2D model Eq. 0 with /i = 0 as a function 
of s — d coupling Bq and a = aa for (1,0) spiral. For com¬ 
parison, the nontrivial phase of ID model with fi — —0.2 and 
—0.4 are shown as the region inside light blue dots and white 
dots, respectively, (b) Energy levels of a system with size 
L = 80a subject to a spiral of wave length A = 27r/|Q| = 16a 
at Bq — 0.3. The orange regions are the qy^s at which the 
topological criterion Eq. Q is satisfied, where one also sees 
the Majorana zero energy edge state, (c) the LDOS along x 
for different values of A/L. Red, green, and blue lines are the 
first, second, and third site away from the edge. The A/L = 2 
case represents a Neel or Bloch domain wall in a ferromagnetic 
insulator. Other parameters are t = —1, A = 0.05. 


The weight of ZBP decreases as the spiral wave length 
A = 27r/|Q| increases. The A/L = 2 case represents a 
Neel or Bloch magnetic domain wall joining two regions 
of opposite spin orientations in a ferromagnetic insulator, 
since it can be viewed as a spiral with half a wave length, 
and the interface to a SC can be described by Eqs. 0 to 
In such case, the ZBP is small but still discernable. 

For a thin film with finite thickness but the s — d cou¬ 
pling only at the interface atomic layer, the Majorana 
state extends over few layers away from the interface, so 
one may need a SC film of few atomic layers thickness 
to observe the Majoranas by any surface probe such as 
scanning tunneling microscope (STM). Even if the mul- 
tiferroic contains domains of different spiral chirality, or 
the spin texture is not perfectly periodic, the Majorana 
edge state still exits at the edge and the boundary be¬ 
tween domains. The large single domain of multiferroics, 
currently of mm size [38], may help to separate the Ma¬ 
jorana fermions over a distance of macroscopic scale. 

SC/skyrmion interface.- Skyrmion spin textures have 
been observed in thin film insulating multiferroics [391 [40] 


FIG. 2: (color online) (a) The spin texture of the pedagogical 
model, which is closely related to a single nonchiral skyrmion 
in (b) a polar lattice and (c) a square lattice, (d) The dis¬ 
persion of the pedagogical model with OBG in f and PBG in 
cp, with Nr = = 40, and (e) if the emerging EM field is 

manually turned off, at Bq = 0.24. The spontaneous current 
pattern at the interface to a single nonchiral skyrmion of size 
9 X 9 is shown in (f), (g) for the normal state (A = 0) and 
(h), (i) for the SG state. Parameters are t = — 1, A = 0.2, 
y = — 0 . 1 . 


at temperatures approaching the typical SC transition 
temperature, with a small magnetic field that presum¬ 
ably has negligible effect on the SC. To gain more un¬ 
derstanding about the SC/skyrmion insulator interface, 
disregarding external magnetic fields, we first consider a 
closely related pedagogical model defined on a square lat¬ 
tice, whose low-energy sector can be studied analytically. 
The spin texture of this model is that shown in Fig.j^a), 
yielding a magnetic field on its interface to an SC 


'Bi = B (sin Oi cos ipi , sin Oi sin ipi , cos Oi ) (9) 


at position (vi^ipi) on the square lattice, where Oi = 
'KTijR and R is the width in r direction. We assume OBC 
along the r direction and PBC along the dp direction. The 
Flamiltonian is described by Eq. M with S = {r, p}. To 
align the spin texture along , one performs the rotation 
in Eq. Q with Ui defined by 


Ui = 


cos ^ — sin^e 

sin^e*"^* cos ^ 


( 10 ) 
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which yields Eq. Q with 

ai^ = 1 + sin^ ^ ^ ’ 

i sin, 

ip 

Pir = e'“^' sin ^ = e'“^'Pir , ( 11 ) 

where Nr and N^ are the number of sites in each direc¬ 
tion. After gauging away the extra phase pis Pis by 
e-t^il 2 g,^ ^icpil 2 g,^ Hamiltomau is 

translationally invariant along (f but not along r because 
Piip (X sin^^ = sm{7rri/R). In the 5 ^ |/i| ^ {t< 5 ,Ao} 
limit, using Eq. (§ of the low-energy effective theory, the 
induced gap along r and cp are proportional to and 
;d*^, and therefore of {pr + ipc^)-wave-like symmetry. 

The spin-conserved aistis and spin-flip Pistis hopping 
in the pia basis contain an emergent electromagnetic 
(EM) field (dTJ [42] coming from the spatial dependence of 
the unitary transformation [43] . This becomes evident in 
the continuous limit Oi 0^ pi p, Ui U ^ and using 

drdrO = 21^jNr <C 1 and d(pdp>p = 27rfN^ ^ 1. 
The {as^Ps} contain the phase gained over one lattice 
constant 

( 7s af' ) = 

where As = i (H/q) WdsU . The factor of i difference 
between Ar and A^ eventually leads to the induced 
(p^+ipc^)-wave-like gap. Eigurej^d) shows the dispersion 
of this pedagogical model, and Eig. [^e) shows the dis¬ 
persion when the emergent EM field is manually turned 
off by setting ai^p = 1 in Eq. ( [TT] ). Without the emer¬ 
gent EM field, the dispersive edge bands expected for the 
induced {pr + ipc^)-wave-like gap are evident, whereas in 
the presence of it the bulk gap is diminished, although 
the trace of edge bands can still be seen in som cases 
(compare Qp ^ 0 and E{n,qp) ~ 0 regions in Eig. j^d) 
and (e)). 

The relevance of this pedagogical model is made clear 
by shrinking the spins at the = 0 edge (green arrows in 
Fig-H a)) into one single spin, which results in a nonchi- 
ral skyrmion on a polar lattice shown in Eig. [^b), with 
the same interface magnetic field described by Eq. 
Certainly these two lattices cannot be mapped to each 
other exactly, but their low energy sectors display simi¬ 
lar features. 

The pedagogical model indicates that the 
SC/skyrmion interface hosts a complex interplay 
between (i) the s-wave gap, (ii) the induced {pr + ipip)- 
wave-like gap, (iii) the emergent EM field, and (iv) the 
s — d coupling Bq. Motivated by the STM generated 
single skyrmion [44] (although with an external magnetic 
field), we proceed to study the SC/skyrmion interface 
on a single open square (SC/SK), whose spin texture 


at position = {xi^yp = |r^| (cos sin(/:?^), as 
shown in Eig.j^c), is that described by Eq. (|^, but with 
Oi = 7 r|r^|/i?(r^). Here, R{Yi) is the length of the straight 
line that passes through connecting the center of the 
square with the edge. The spontaneous current at site i 
can be calculated from Eq. 0 by 

J* = ■ (13) 

(tS 

As shown in Eig. Hf) and (g), in the normal state 
(NM/SK) there is a persistent current whose vorticity 
strongly depends on the emergent EM field, the s — d 
coupling Ho, and finite chemical potential. Because the 
persistent current also has edge component, the topolog¬ 
ical edge current alone is hard to be quantified in the 
SC state from the pattern shown in Eig. [^h) and (i). 
At Ho ^ A, the current pattern in the SC state recov¬ 
ers that of the NM state, as can be seen by comparing 
Eig. I^g) and (i). Eor the chiral skyrmion seen in most 
experiments, the form of H^ and H^ in Eq. 0 are ex¬ 
changed, which gives the same result as and can be 
trivially exchanged for spin-independent quantities such 
as J^. These results suggest a vortex-like state at the 
NM/skyrmion lattice or SC/skyrmion lattice interface, 
whose vorticity depends on material properties. 

In summary, we propose that a broad class of non- 
collinear magnetic orders, including a large part of those 
discovered in multiferroic insulators, and the Bloch and 
Neel domain walls in collinear magnetic insulators, can 
be used to practically generate Majorana edge states at 
their interface to a conventional SC. The advantages of 
these systems include a much larger parameter space to 
stabilize the edge states compared to ID proposals, the 
longer range magnetic order may help to separate the 
edge states over a macroscopic distance, and adjusting 
chemical potential is not necessary. The proximity to a 
skyrmion induces an inhomogeneous {pr + ipc^)-wave-like 
pairing in the SC under the influence of an emergent elec¬ 
tromagnetic field, and consequently a vortex-like state 
that features both a bulk persistent current and an edge 
current. 

We thank P. W. Brouwer, Y.-H. Liu, and E. von Oppen 
for stimulating discussions. 
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